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Abstract 

We show that there is a query expressible in first-order logic over the 
reals that returns, on any given semi-algebraic set A, for every point a 
radius around which A is conical. We obtain this result by combining 
famous results from calculus and real algebraic geometry, notably Sard's 
theorem and Thorn's first isotopy lemma, with recent algorithmic results 
by Rannou. 



1 Introduction 

The framework of constraint databases, introduced by Kanellakis, Kuper and 
Revesz || , provides a nice theoretical model for spatial databases fTTJ . A spatial 



dataset is modeled using real polynomial inequality constraints; such sets are also 
known as semi-algebraic sets |], |3|] . The relational calculus (first-order logic) with 
real polynomial constraints then serves as a basic query language, denoted here 
as FO. 

The study of the expressive power of query languages for constraint databases 
is an active domain of research | XOf] . One of the problems in particular that 



received attention in recent years is that of determining the exact power of FO 
in expressing topological properties of spatial databases P, |7|, p|, |I3[1 . One such 



property, which is central in this research, is that locally around each point, a 
semi-algebraic set has the topology of a cone. A radius at which this behavior 
shows is called a cone radius around the point for the set. 

Accordingly, a cone radius query is a query that returns, for a semi-algebraic 
set A, a set of pairs (p, r) giving for every point pa. cone radius r in p for A. In this 
paper, we show that there exists an FO formula expressing a cone radius query. So 



far, this was only known in two dimensions |J. Expressibility of the cone radius, 
apart from being a natural question in itself, has also applications: for example, 
it has been linked to the expressibility of piecewise linear approximations j|. 



2 Preliminaries 

2.1 Spatial databases and Queries 

A semi-algebraic set in W n is a finite union of sets definable by conditions of the 
form 

fi(x) = f 2 (x) = ■■■ = f k (x) = 0, gi(x) > 0, g 2 (x) > 0, . . . , g e (x) > 0, 

with x = (xi, . . . ,x n ) G E n , and where f\ (x),... , fk(x), gi(x), . . . , gi(x) are 
multi-variate polynomials in the variables with real coefficients. A 

database schema S is a finite set of relation names, each with a given arity. A 
database over S assigns to each S G S a semi-algebraic set S D in MP, if n is the 
arity of S. A k-ary query over S is a function mapping each database over S to 
a semi-algebraic set in M. k . 

As query language we use first-order logic (FO) over the vocabulary (+, - ,0,1, <) 
expanded with the relation names in S. A formula tp(xi, . . . ,Xk) expresses the 
fc-ary query defined by 

tp(D) := {(at, . . . , a k ) e R k | (R, D) |= ip( ai , ... ,a k )}, 

for any database D. Note that <p(D) is always semi-algebraic because all rela- 
tions in D are; indeed, by Tarski's theorem the relations that are first-order 
definable on the real ordered field are precisely the semi-algebraic sets. 

An example of a query expressed in FO is the following: Let S be a schema 
containing the relation name S. Consider the FO-formula 

(Pmt(x) ■= (3e > 0)(Wx[) ■ ■ ■ (Wx' n )(\\x - x\\ < e -> S(x[, . . . ,x' n )). 

For any database D, tp- mt (D) equals the interior of S D . However, not every 
query is first-order expressible: the query which asks whether a set is connected 
is not expressible in FO. This result and other results related to constraint 



databases have recently been collected in a single volume 1C . 



2.2 Cones 

Let A C W 1 be a semi-algebraic set and p e R n a point not in A. We define 
the cone with base A and top p as the union of all closed line segments between 
p and points in A. We denote this set with Cone(A,p) := {tb + (1 — t)p \ b G 
A, ^ t ^ 1}. For a point pGl" and e > 0, denote the closed ball centered at 
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p with radius e by B n (p,e), and denote the sphere centered at p with radius e, 
by S 11 ^ 1 ^, e). The following well-known theorem says that, locally around each 
point, a semi-algebraic set has the topology of a cone. 

Theorem ([l], |3[]). Let A C R n be a semi- algebraic set and p a point of A. 
Then there is a real number e > such that the intersection A D B n (p,e) is 
homeomorphic to the set Cone(^4 fl S n ~ 1 (p,e),p). 

Any real number e > as in the lemma is called a cone radius of A in p. 

Let 5 be a schema containing a relation name 5* of arity n. A cone radius 
query Qradius is a query which maps any database D over S to a set of pairs 
(p, r) G R n x R such that for every point p G S D there exists at least one pair 
(p,r) G Qradius(-D), and for every (p,r) G Qradius(-D), r is a cone radius in p for 

We will use the following notation: Let A C S C R n , the closure of A in 
B is denoted by c1b(^4), and ints(A) indicates the interior of A in 5. When 
the ambient space B is R n , we omit the subscript B. We denote cl(A) — A (the 
frontier of A) with OA. 



2.3 Whitney stratification 

Every semi-algebraic set ACR" can be "nicely" decomposed in a finite sequence 
Z of n + 1 semi-algebraic sets Z , . . . ,Z n , called strata, with the following prop- 
erties. For each i = 0, . . . , n: 

1. Zi is either a C 1 semi- algebraic set in R n of dimension i, or an empty set; 
and 

2. each triple (Zi,p, Zj) for i < j and p & Zi has the Whitney property. 

Such a decomposition is called a C 1 - Whitney stratification of A. We will not need 
the precise definitions of when a set is C 1 and of the Whitney property. We refer 



to the paper of Rannou |12| for more details. 

We remark that in this paper we do not require the frontier condition, which 
says that the frontier of a stratum is the union of lower dimensional strata, and 
also do not suppose strata to be connected. Both properties are not necessary for 
Thorn's first isotopy lemma R [TJI], which we will use in our proof in Section 01. 



3 Constructing a Whitney stratification 

Let A be a semi-algebraic set in R n . We shall construct a C 1 - Whitney stratifica- 
tion Z of the closure c\{A) such that A is the union of connected components of 
strata of Z. We then say that Z is compatible with A. This construction will be 
expressible in FO. 
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Our construction is an adaptation of the construction given by Shiota fll4 , 
Lemma 1.2.2]. 

We define Z n as the subset of A where A is locally C 1 and of dimension n. 
Now suppose that the strata Z n , . . . , Z k +i have already been constructed. Then 
the stratum Z k is constructed as follows. Define A = A and Ai = OA. For 
i = 0, 1 construct 



R % k := {p G Ai — Zj | Ai is C 1 and of dimension k 



j=k+l 



(1) 



in a neighborhood of p} 



W l k := p| int^({p e R\ | (R\,p, Zj) has the Whitney property}) (2) 

j=k+l 

Zl:=Wl-d{R\-*). (3) 

Then we define Z k : = Z'u Z\. 

The stratum Z k has indeed the desired properties: By (Jlj) it is C 1 and of 
dimension k, (0) guarantees that for all points in Z k , and for any j > k, the 
triples (Z k ,p, Zj) have the Whitney property, and (H) ensures that the connected 
components of Z k lie either completely in A or completely in dA. 



It is well known [fT7|, [LJ] that the dimension of Aj — IJj=fc %j 1S strictly less 
than the dimension of Ai — U"_ fe+1 Zj for i = 0,1. Hence, the stratification Z 
will consists of exactly n + 1 strata Zk, some of which may be empty. 

We now show that the above construction is in FO. 

Theorem 1. Let S be a database schema containing a relation name S of arity 
n. The n-ary query Qk-stratum, which takes as input a database D over S, and 
returns the kth stratum Z k of the stratification Z constructed above for A = S D , 
is expressible in FO. 

In order to prove FO-expressibility of these queries, it is sufficient to show 
that the sets R l k , W k , and Z\ occurring in the construction are FO-expressible. 
But this follows immediately from the work of Rannou [|T2"|| . Indeed, from that 
work we can deduce the following lemma, which immediately implies Theorem |l} 

Lemma 1. (i) Let S be a database schema containing a relation name S of arity 
n. The n + 1 queries of arity n, defined as 

Qk-reg(D) := {f G R" | i G S D A (S D is C 1 and of dimension k 

in an open neighborhood of x)}, 

for any database D over S, are all expressible in FO. 

(ii) Let S be a database schema containing two relation names Si and S2 of arity 
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n. The n-ary query, defined as 

Qwhtiney(D) := {x G R n | Sf, are C 1 and (Sf , f, S 2 D ) 

aas t/ie Whitney property}, 

for any database D over S, is expressible in FO. 

4 Expressing the cone radius in FO 

We are now ready to prove the main result of this paper. 

Theorem 2. There exists an FO -expressible cone radius query. 

Proof. Consider a semi-algebraic set A in R n , and let Z be the C 1 - Whitney 
stratification of cl(A) compatible with A. Let p G cl(A) and define the C 1 -map 

fy: c\(A) -»R:ih ||f -p1| 2 . 



We will need Thorn's First Isotopy Lemma [13]. Applied to the C 1 -map and 



the C 1 -Whitney stratification 2, this lemma can be stated as follows: For any 

a < 6, 

(a) If fp is proper, i.e., f# ([a, b]) is compact, and 

(b) if for each stratum Z £ Z, the restriction 

fy\ (Z H int(B n (p, b) - B n (p, a))) -> (a, 6) CR 

has no critical points (this will be explained later), 
then for any c G (a, 6), there exists a homeomorphism 

/i : cl(A) n int(B n (p, b) - B n (p, a)) -> (cl(A) n S^p, c)) x (a, 6). 
Moreover, this homeomorphism satisfies the following two properties: 

(i) For each d G (a, 6), h(d(A) n S ,n - 1 (p, d)) = (cl(A) n S^ 1 ^ c)) x {d}, and 

(ii) /i(Znint(i3 n (p, b)—B n (p, a))) = (ZnS' n_1 (p, c)) x (a, 6) is a homeomorphism 
for every connected component Z of a stratum of i?. 

This statement of Thorn's First Isotopy Lemma is a specialized form of Theorem 



II.6.2 in Shiota [|TJ. 

Remark that condition (a) is automatically satisfied. Indeed, the inverse 
image by f$ of any interval [a, b] C R is equal to cl(A) fl (B n (p, b) — int (B n (p, a))) 
which is closed and bounded in R". 
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Claim 1. If condition (b) is satisfied for < b (so a = 0), then every c G (0, b) 
is a cone radius of A in p. 

Proof of Claim. Take an arbitrary real number c G (0,6). The lemma gives a 
homeomorphism 

h : cl{A) n intC^Cp, b) - {p}) -> (cl(A) n S n ~\p, c)) x (0, 6). 

By property (i), we obtain a homeomorphism 

/i! : d(A) n (iT(p, c) - {p}) - (cl(A) n S n ~\p, c)) x (0, c], 

which equals the restriction /i |cl(A) D (B n (p,c) — {p}). Since the cylinder 
(cl(A) n S n ~ l (p, c)) x (0,c] is clearly homeomorphic to the cone Cone(cl(A) fl 
S n ~ 1 (p,c),p) — {p}, e.g., by the homeomorphism 

h 2 (x,t) := (1 )p + (-)x forte(0,c], 

c c 

we obtain a homeomorphism 

ha := /i2 o hi : d(A) n (5 n (p, c) - {p}) -> Cone(cl(A) n S n -\p, c),p) - {p}. 

It is easily verified that h 2 : (ZnS™- 1 ^ c)) x (0, c] -> Cone(Zn5 n ~ 1 (p, c),p)- 
{p} is a homeomorphism for each connected component Z of a stratum of i?. 
Since /ii also satisfies property (ii), we have that h% : (Z n (B n (p,c) — {p}) — > 
Cone(Zn S n ~ 1 (p, c),p) — {p} is a homeomorphism for each connected component 
Z of a stratum of J?. 

This implies that the restriction /z = /i 3 |yl fl (B n (p, c) — {p}) is a homeomor- 
phism from /i(^4 fl (B n (p, c) — {p}) to Cone(A fl S^ 1 ^, c),p) — {p\, because A is 
the union of connected components of strata of Z. 

The homeomorphism h can easily be extended to the point p, hence c is indeed 
a cone radius as desired. □ 

Let 5 be a schema containing a relation name S of arity n, and let D be a 
database over S. By Claim [l], we can define the following cone radius query 

QradiusP) := {(P, r) E W l x M \ p E S D and r G (0, 6)}, 

where the interval (0, b) satisfies condition (b) for the map fp and semi-algebraic 
set A = S D . Let us express this query in FO. 
We define the critical point query as 

Q crit (D) := {(p, x) G R n x R n | p G S D and x G Qfc- S tratum(£>) for a certain 

and of is a critical point of fp restricted to Qk-stTa,tum(D)}. 

The critical points of restricted to a stratum Z, are the points feZ where 
the differential map dg{fp\Z) is not surjective. 
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Claim 2. A point x 6 I" is a critical point of fp\Z if and only if the tangent 
space of Z in x is orthogonal to the vector x — p. 

Proof of Claim. We compute the differential d^{f^\Z) as follows: Locally around 
x, we may assume that the projection on the first k coordinates II : Z — > U C 
IR fc , is a homeomorphism, where k is the dimension of Z in x. By definition 
of the differential, d s (f$\Z) = (d {xu ... !Xk) g)(d (xu ... ^U' 1 )- 1 , where g = (f\Z) o 
II -1 . By the C 1 Inverse Function Theorem, we may assume that II -1 : U — ► 
Z : (xi, . . . , Xk) >->• {xi, ... ,x k , cpk+i, ... , ip n ), where <#(xi, ... ,x k ) are C 1 -maps, 
and hence g : U \->R: (xi, ... ,x k ) — J2i=i( x i ~ Pi) 2 + E^fc+i^jO^i) • • • > x k) ~ 
Pj) 2 - An elementary calculation shows that the differential of f^\Z in x is the 
vector 

d x (fp\Z) = 2 || (xi -pi) + E^-Pi)^^-^^) » °»- >° 

\\ i=fc+l ' i=l,...,k n — k times/ 

Since gL^...^)!! -1 is an isomorphism between the tangent space T( Xlr .. <Xk )U of 
in the projection H(x), and the tangent space T$Z of Z in x, any tangent vector 
(v±, . . . ,v n ) G TgZ is of the form (d^x 1> ... >Xh )H~' L )(vi, . . . ,v k ). More specifically, 
any tangent vector v G T$Z can be written as 

k k 

(«!,... ,v„) = (ui,... g fc+1 (^i,--- ,x fc )ui,... ,x fc )uj). 

i=i 1 i=i % 

Hence, the product 



• . . - j X k IVi I , 



fe n / k Q 

dx(fp\z) -v = 2 ^2(xi - p^ + 2 ( x j - Pi) I E 

i=l j"=Jfc+l \i=l * 

is equal to 2Y2i=i( x i ~ Pi) v i- This implies that the differential map dg(f^\Z) is 
not surjective if and only if 2 Yli=i( x i ~Pi) v i = f° r & U tangent vectors u G T^Z. 
This proves the Claim. □ 

The proof of the theorem now continues as follows. The tangent space query 
Qtangent(^) := {(x, ??) 6 1" x 1" S D is C\ x G 5° and V G T^}, 
is expressible in FO |TJ , Lemma 2] . Because the orthogonality of two vectors can 



be easily expressed in FO, the formula 

n 

<A:rit(P, X*) := Stf) A \/ W (^tangent (^-stratum (5') )(X, ->• (x - p) • U = 0) 
j=0 
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expresses Q cr i t correctly by Claim 2. Here, ^-stratum denotes an FO-formula 
expressing Qj-stratum for j = 0, . . . ,n, and (/'tangent is an FO formula expressing 

Qtangent ■ 

A critical value of fp is the image by fp of a critical point. The query which 
returns the set of critical values is expressible in FO by the formula 

tp v3 .\(p, r) := 3x (y?crit(p, x) Ar = fp(x)) . 

We now observe that {r £ R | (p va \(p,r)} is finite for each p. Indeed, the set 
of critical points {x £ M. n \ ip CTit (p,x)} is semi- algebraic and hence admits a C 1 - 
cell decomposition C = {Ci, . . . ,C m } such that /|Cj is C 1 Jl6|]. Sard's Theorem 
for C 1 -maps |18| implies that each fp\Ci attains only a finite number of values. 
Hence the image by fp of the set of critical points is finite. 

This implies that either there are no critical values, or there exists a minimal 
critical value. By Claim 1, any value smaller than this minimal value is a cone 
radius. We therefore conclude that the query expressed in FO as 

^radius (p,r) := (Vr')(v?vai(p>0 -> r < r'), 

is a cone radius query, as desired. □ 

Acknowledgement 

The author wishes to thank Jan Van den Bussche for many interesting discussions 
and critical remarks, which led to improvements in the presentation of the paper. 



References 

[1] R. Benedetti and J.J. Risler. Real Algebraic and Semi- algebraic Sets. Her- 
mann, Paris, 1990. 

[2] M. Benedikt, M. Grohe, L. Libkin, and L. Segoufin. Reachability and con- 
nectivity queries in constraint databases. In Proceedings of the 19th ACM 
Symposium on Principles of Database Systems, pages 104-115. ACM Press, 
2000. 

[3] J. Bochnak, M. Coste, and M.-F. Roy. Real Algebraic Geometry, volume 36 
of Ergebenisse der Mathematik und ihrer Grenzgebiete. Springer, 1998. 

[4] F. Geerts and B. Kuijpers. Linear approximation of planar spatial databases 
using transitive-closure logic. In Proceedings of the Nineteenth ACM Sympo- 
sium on Principles of Database Systems, pages 126-135. ACM Press, 2000. 



8 



[5] F. Geerts and B. Kuijpers. Expressing topological connectivity of spatial 
databases. In A. Mendelzon and R. Connor, editors, Research Issues in 
Structured and Semistructured Database Programming, volume 1949 of Lec- 
ture Notes in Computer Science, pages 224-238. Springer- Verlag, 2001. 

[6] C.G. Gibson, K. Wirthmiiller, A. A. du Plessis, and E.J.N. Looijenga. Topo- 
logical Stability of Smooth Mappings, volume 552 of Lecture Notes in Math- 
ematics. Springer, 1976. 

[7] M. Grohe and L. Segoufm. On first-order topological queries. In Proceedings 
of the 15th IEEE Symposium on Logic in Computer Science, pages 349-360. 
IEEE Press, 2000. 

[8] P.C. Kanellakis, G.M. Kuper, and P.Z. Revesz. Constraint query languages. 
Journal of Computer and System Sciences, 51(l):26-52, 1995. 

[9] B. Kuijpers, J. Paredaens, and J. Van den Bussche. Topological elemen- 
tary equivalence of closed semi-algebraic sets in the real plane. Journal of 
Symbolic Logic, 65(4):1530-1555, 2000. 

[10] G.M. Kuper, J. Paredaens, and L. Libkin, editors. Constraint Databases. 
Springer- Verlag, 2000. 

[11] J. Paredaens, J. Van den Bussche, and D. Van Gucht. Towards a theory 
of spatial database queries. In Proceedings of the 13th ACM Symposium on 
Principles of Database Systems, pages 279-288. ACM Press, 1994. 

[12] E. Rannou. The complexity of stratification computation. Discrete and 
Computational Geometry, 19:47-79, 1998. 

[13] L. Segoufin and V. Vianu. Querying spatial databases via topological invari- 
ants. Journal of Computer and System Sciences, 61(2):276-301, 2000. 

[14] M. Shiota. Geometry of Subanalytic and Semialgebraic Sets. Birkhauser, 
1997. 

[15] L. van den Dries. Alfred Tarski's elimination theory for real closed fields. 
Journal of Symbolic Logic, 53(1):7-19, 1988. 

[16] L. van den Dries. Tame Topology and o-minimal Structures. Cambridge 
University Press, 1998. 

[17] L. van den Dries and C. Miller. Geometric categories and o-minimal struc- 
tures. Duke Mathematical Journal, 82(2):497-540, 1996. 

[18] A.J. Wilkie. A theorem of the complement and some new o-minimal struc- 
tures. Selecta Mathematica, New Series, 5:397-421, 1999. 



9 



